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1. Motivation 

The Monte Carlo simulation of lattice field theories at finite chemical potential /i is a notoriously 
difficult task due to the complex action problem: At nonzero /r the action S is complex and the 
Boltzmann factor cannot be used as a probability weight. A powerful approach to overcome 
the complex action problem is to find a re-paramefrizafion of fhe parfifion sum Z in ferms of new 
variables, such fhaf Z is a sum over only real and posifive confribufions. Then a Monfe Carlo 
simulafion can be sef up in terms of fhese so-called dual variables, which are loops for maffer fields 
and surfaces for gauge fields. The dual approach has been successfully implemenfed for several 
inferesling lattice field fheories af finite densify (see, e.g., [1, 2, 3] for examples). 

So far fhe dualizafion was mainly resfricfed fo bosonic fheories, since fhe fermions con- 
fribufe addifional minus signs for fhe maffer loops, which come from fhe Pauli principle and fhe 
fraces over fhe y-mafrices. However, infegrafing ouf fhe gauge fields leads fo addifional consfrainfs 
for admissible dual configurations of loops and surfaces, which can lead fo cancellafions of fhe 
fermion signs. In order fo sfudy fhese emerging dual sfrucfures we here analyze lattice QED wifh 
fwo flavors of Wilson fermions and a chemical pofenfial. For fhis model fhe gauge fields can be 
dualized complefely (fhis is possible for all abelian gauge fheories) and one can focus on sfudying 
fhe sfrucfure of fhe fermion loops inferacfing wifh fhe gauge surfaces. In fhis confribufion we sfudy 
fhe loops as obfained from a partially summed hopping expansion and fheir inferacfion wifh fhe 
dual gauge surfaces. We show fhaf for cerfain classes of loops (quasi-planar loops), fhe fofal signs 
are posifive and we presenf resulfs for a dual simulafion af finite chemical pofenfial. 


2. Definition of the model and hopping expansion of the fermion determinant 

In fhe conventional represenfafion fhe partition sum of lattice QED wifh fwo flavors is given by 

Z = j 9[U] detD[U] detD[U*] . (2.1) 

Here fhe fermions are already infegrafed ouf and only fhe pafh integral over fhe link variables 
Uv{n) G f/(l) remains, where &[U] denotes fhe producf of fhe Haar measures for all gauge links. 
For fhe gauge acfion we use fhe Wilson form, = — j8 Y,p Re Up, where j3 denofes fhe inverse 

gauge coupling, fhe sum runs over all plaqueffes p, and Up is fhe producf of link variables along 
fhe confour of p. We consider fwo mass-degenerate flavors of opposife charge, such fhaf Gauss’ 
law is obeyed. As a consequence fhe conjugafe gauge links Uy{n) are used in fhe second fermion 
determinant. The Wilson Dirac operator is denoted by D[f7] = 1 — kH[U], with the standard Wil¬ 
son hopping matrix H[U]{n,m) = Y,v Uv{n) 5„+v,m + Uy{n — v)* > 

where K is the hopping parameter and /t the chemical potential. Both parameters are chosen to have 
the same values for the two flavors (but remember that for the second flavor the complex conjugate 
link variables are used in the hopping matrix), n and m label the sites of a A| x Nj lattice with 
boundary conditions that are periodic in the three spatial directions and anti-periodic in time. 

In order to obtain the fermion loops of the dual representation we apply hopping expansion 
of the determinant (see, e.g., [4] for an elementary introduction), i.e., an expansion in small K 
(large mass). Using the trace-log formula the fermion determinant can be rewritten as detD = 
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exp(—^ Tr //”), where we have already expanded in K the logarithm in the exponent. The 
traces Tr H” over powers of the hopping matrix give rise to non-backtracking closed loops of 
length n which are dressed with the corresponding hopping terms along their contour. Thus the 
^ Tr can be viewed as a sum over all closed loops. Among the set of all loops there are also 
those where a simple loop is iterated several times. Thus we can split the sum into a sum over all 
non-iterating loops and a sum over the iterations of the loops. The sum over the non-iterating loops 
in the exponent we write as a product over the exponentials, and we keep in the exponent only the 
sum over the iterations. However, this sum over iterations in the exponent can be written again as 
a determinant over only the Dirac indices of the product of hopping terms along the contour of the 
loop. Thus we have brought the determinant into the form 

detD= det(l-c/r/)= (l-tr(r/)c/-hsumdet(r;)c^) . (2.2) 

Here .ifNi is the set of all oriented, non-backtracking, non-iterating loops /, T/ is the product of the 
matrices (1 =F 7v)/2 along the loop I, and c/ is given by 

Cl = _ (-2.3) 

By 1T(/) we denote the number of times the loop I winds around compact time (here we obtain an 
extra minus sign from every winding due to the anti-periodic temporal boundary conditions of the 
fermions), and Ui is the product of the link variables along the loop /. 

The determinant in (2.2) now has the form of a product over loops in where each loop 
/ can contribute in three ways: The first possibility is a factor of 1, which is what we refer to as ’’not 
activated”. However, the loop I can also contribute a factor of —tr(r/)c/ ("single occupancy"), or 
finally, / confribufes a factor sumdef(r/) c^, which we refer to as "double occupancy" of fhe loop 1 . 

In fhe second sfep of (2.2) we have wriffen fhe deferminanfs for fhe individual loops I wifh 
fhe characferisfic polynomial ;^c/r/(■^) = del( A—c/T;) af A = 1. A characferisfic polynomial can be 
expanded in a sum over deferminanfs of reduced mafrices (some columns and rows deleted). Since 
fhe T/ are producfs of fhe rank-2 projectors (1 =F Jv) /2, only terms where fwo or fhree columns and 
rows are delefed remain in fhaf expansion. This leads fo fhe rhs. of (2.2), where we infroduced 
sumdef (A) = Y,<n.m> delA<„ where fhe sum runs over all ordered pairs <n,m> wifh 1 < n < 
m < 4, and A^n.m> is the reduced matrix where the n-th and m-th rows and columns are removed. 

From (2.2) and the form of the c/ it is obvious that only temporally winding loops couple to 
the chemical potential. Thus we now restrict the sum over loops in (2.2) to winding loops, where 
for now we do not distinguish between winding around compact time or space directions. We 
consider so-called quasi-planar loops, which are loops that after one or more steps in a direction 
perpendicular to the winding direction again continue in the winding direction. A key step is the 
evaluation of the signs coming from the Clifford algebra, which in our presentation corresponds 
to computing tr (T/) and sumdet(r/). It has been known for a long time that for planar loops I the 
traces over the F/ can be computed in closed form [5]. Using similar techniques as in [5] one finds 
for quasi-planar loops 

tr(r/)=(0 and sumdef (F,) = , (2.4) 

where t denofes fhe number of deviafions from fhe main direction. 


3 



Dual simulation of finite density lattice QED at large mass 


Christof Gattringer 


3. Integrating out the gauge fields and final form of the dual representation 


The second step in the dualization is to integrate out the gauge fields. They come from two sources: 
The fermion loops are dressed with gauge link variables along their contour and in the plaquettes 
of the gauge action we have products of four link variables along all unit squares of the lattice. The 
important key formula for integrating the link variables is (the jv{n) are arbitrary integers) 


n^vW 


jv{n) 


ny 


ny 


(3.1) 


i.e., only contributions where the links along the fermion loops are compensated with links from 
the plaquettes of the gauge action do survive (see also the discussion in [3]). To organize the gauge 
integration we expand the Boltzmann factor of the gauge action using the well known series 

oo 

,-som = (>.',+1';) = pj ^ u "/, (3.2) 

P P 


where /|„^| (jS/2 ) denotes the modified Bessel funcfions, and the integers rip are referred to as the 
plaquette occupation numbers for the plaquette p. These plaquette occupation numbers describe 
the dual degrees of the U(l) gauge fields. The series in (3.2) has ferms wifh arbifrary infeger powers 
lip of a plaqueffe Up, where negative powers correspond fo an orienfation of the plaquette in the 
mathematically negative sense. Only those configurations of the plaquette occupation numbers iip 
survive, where the occupied plaquettes compensate the links along the fermion loops. This implies, 
that for obtaining non-vanishing dual configurations one builds up surfaces of occupied plaquettes 
p (i.e., rip / 0), and the surfaces are either closed or bounded by the fermion loops. Note that for 
the second flavor the link variables along the loops are complex conjugate and plaquettes with the 
opposite orientation have to be used to compensate that flux. 

With (3.1), (3.2) and the results for the loops from the previous section we obtain 


Z = 


L 

{h,h,n} 


gMW+W,) ^Li+Lj 

2fl+^2 


n'Kl(/5/2)^ 


(3.3) 


The sum is over the sets {fi},{f 2 } of all quasi-planar, non-iterating loops l\ and I 2 representing 
the dual fermion degrees of freedom for the two flavors, where each loop can have either single or 
double occupation. By W \, W 2 we denote the corresponding total winding numbers for all loops of 
the two flavors. L\ and L 2 denote the total length of all the loops of the two flavors, and T\ and T 2 
are the shape factors of the loops, computed from the number of deviations of the loops from the 
winding direction according to (2.4), with the trace being used for single occupation and sumdet for 
double occupation. In Z we also sum over the dual degrees of freedom for the gauge fields: We sum 
over all admissible configurations {n} of the plaquette occupation number iip, where ’’admissible” 
refers to the fact that the total flux along the fermion loops is compensated by plaquettes. As 
mentioned before, this amounts to summing over surfaces made from occupied plaquettes, where 
the surfaces are either closed or bounded by fermion loops (see also the discussion in [3]). 

An important constraint comes from Gauss’ law which requires an equal number of loops 
that have their links oriented forward and backward in time (double occupied loops count twice). 
In the dual language this is manifest in the fact that only for vanishing link flux from the loops we 
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can obey all constraints with finitely many occupied plaquettes. These constraints can be obeyed 
for example by a forward and a backward winding fermion loop of the same flavor, or by two loops 
of different flavor but with the same orientation. In all cases the total number of occupied loops 
is even (double occupation counts twice) and we do not have any remaining minus signs in (3.3) 
from single occupation of loops or the anti-periodic temporal boundary conditions of the fermions. 

4. Dual simulation of the model and numerical results 

We consider the model described by (3.3) in two variants: The first variant uses the full dynamics 
of the quasi-planar loops, while in the second variant we restrict the dual fermionic degrees further 
and only consider linear (straight) loops winding in time direction. The latter case corresponds to 
static charges and we refer to it as the "linear loop model". 

Both variants of the model are updated with local Monte Carlo steps. For the dual gauge 
degrees of freedom we use two steps that can update also pure gauge theory: 1) We propose a 
change of all plaquette occupation numbers Up of an elementary 3-cube by ±1, arranged such 
that the net flux on all edges of the cube remains unchanged. This proposal is accepted with a 
Metropolis step according to the weights in (3.3) and leads to an admissible configuration. 2) We 
propose a change of all plaquette occupation numbers rip on one of the 6 coordinate planes by ± 1. 
This again leads to an admissible configuration which we accept with a Metropolis step. 

For the linear loop model two updates of the fermion loops are used: 1) We insert pairs 
of temporally winding linear loops at the same spatial position, such that the constraints are kept 
intact. This can either be a pair of oppositely oriented loops of the same flavor, or two loops with 
the same orientation but with different flavor. In the Metropolis acceptance step of the latter, also 
the chemical potential enters. 2) We allow a linear loop to hop in space. In that case we need to 
change the plaquette occupation numbers rip in a strip that connects the old and the new position of 
the loop, and the corresponding weights determine the Metropolis acceptance. Finally, for updating 
the full quasi-planar model, we also consider local, quasi-planar deformations of the loop. Here 
the Metropolis acceptance is computed from the changing length and shape of the loop, as well as 
from the weights of the plaquette occupation numbers that have to be changed in order to remain 
in the set of admissible configurations. 

Various bulk observables were considered in our study - they can be simply obtained as 
derivatives with respect to the parameters jS, K and /r applied to the conventional and the dual 
representations. Here we restrict ourselves to presenting results for the particle number density 
« = 1 /N^Nrd InZjdn, which we use for studying condensation phenomena. 

The simulations were performed on various lattice sizes with Ns and Nt ranging from 4 to 
16. We typically used 0(10'*) sweeps combining all updates discussed above for equilibration and 
the configurations used for measuring the observables are separated by 0(10) combined sweeps 
for decorrelation. The statistics for the observables is typically 0(10^*) measurements and we 
show statistical errors evaluated with the jackknife method. For the results presented here we used 
K = 0.2 and compared the behavior at j3 = 0.9 (confined phase) with j3 = 1.1 (Coulomb phase). 

In Fig. 1 we show the results for n as a function of /i for both the linear loop model (top 
row of plots) and the quasi-planar model (bottom). We compare the results for different spatial 
volumes to study the emergence of a condensation transition. We observe such a transition in both 
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Linear loops, p = 0.9, k = 0.2 



Linear loops, p = 1.1, k = 0.2 



Qnasiplanar loops, p = 0.9, k = 0.2 Quasiplanar loops, p = 1.1, k = 0.2 




Figure 1: The fermion number density n as a function of the chemical potential. The top row is for the linear 
loop model, while the bottom is for quasi-planar loops. On the Ihs. we show the results for p = 0.9, i.e., the 
confining phase, while the rhs. is for the Coulomb phase (j3 = 1.1). The hopping parameter is set to if = 0.2 
and we compare the results for different spatial volumes. 


the confining and the Coulomb phases of both models, although at different values of /i. A striking 
difference is that for the linear loop model we clearly observe saturation of the particle density in 
the condensed phase at a value of n = 2, as is expected since each loop can be occupied at most 
twice. For the quasi-planar loops we do not observe saturation. This is due to the fact that there is an 
infinity of quasi-planar loops that can appear in the expansion of the fermion determinant. On the 
other hand we know that the net particle number is bounded by 2A|, which is the largest power of 
the fugacity that can appear in the fermion determinant. This shows that subtle cancellations of 
terms from a hopping expansion of the fermion determinant are an important aspect that should be 
understood when using this expansion for a finite density quantum field theory. In Fig. 2 we show 
for the linear loop model the results for the density n as a function of /f at a fixed spatial volume but 
now for varying Nj, such that the temperature changes at a fixed physical scale. As expected for 
a condensation transition, we find that for larger temperature (smaller Nj) the transition becomes 
more washed out. A similar behavior was observed for the quasi-planar loops. 

In this exploratory study we analyzed the hopping expansion approach for finding a dual 
description of fermions a finite chemical potential. In particular we were interested in a case where 
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Linear loops, p = 0.9, k = 0.2 Linear loops, p = 1.1, k = 0.2 




Figure 2: The fermion number density « as a function of the chemical potential for the linear loop model. 
We compare the results for different temperatures T, i.e., we vary the temporal extent at a fixed spatial 
volume of 16^. On the Ihs. we show the results for = 0.9, i.e., the confining phase, while the rhs. is for the 
Coulomb phase (j3 = 1.1). The hopping parameter is set to K = 0.2 and we compare different Nj■ 

the fermions can be completely dualized, such that effects of the fermion loops can be studied 
without artifacts from, e.g., a strong coupling expansion. We find that the complete dualization of 
the gauge fields leads to interesting cancellations of minus signs due to Gauss’ law. Furthermore 
we demonstrate that the appearance of saturation in the condensed phase depends on the truncation 
used for the hopping expansion. 
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